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Scales engagement can contribute significantly to nonlinear bending behavior of elastic substrates 
with rigid biomimetic scales. In this letter, we investigate the role of friction in modulating the 
nonlinearity that arises due to self-contact of scales through an analytical investigation. We model 
the friction as dry Coulomb type friction between rigid links and the substrate is taken to be linear 
elastic. Our results reveal that frictional effects give rise to two possible locking mechanisms, namely 
static friction lock and kinetic friction lock. These locks arise due to a combination of interfacial 
behavior and geometry. In addition to these extremes, the frictional behavior is found to increase 
stiffness of the structure. This dual nature of friction which influences both system operation and 
its terminal limit results in the maximum relative frictional work to lie at intermediate friction 
coefficients and not at the extremes of frictional limits. 

PACS numbers: 


Dermal scales can provide considerable mechanical mm 
and multifunctional advantages HD without adding sig¬ 
nificant weight to organisms making them ubiquitous in 
the animal kingdom. These properties along with their 
relative simplicity make them attractive candidates for 
adoption into synthetic materials development paradigms 
in order address the challenges that arise from balanc¬ 
ing weight, mechanical performance and functionality. 
To this end, elementary structures mimicking the nat¬ 
ural scale-systems can be fabricated through partially 
embedding plate like surface appendages in a soft sub¬ 
strate [Ms]. The analysis of these systems which ex¬ 
hibit rich mechanical behavior indicate that the primary 
source of complexity is the engagement of scales with 
each other [ms]- This mechanical enhancement, tai- 
lorable through a straightforward variation of geometry, 
could be leveraged for material systems which require 
tunability of elastic behavior and energetic landscape. 
For a compliant beam like substrate, this scale engage¬ 
ment , which can occur due to change of curvature of the 
parent beam, induces remarkable nonlinearity in bending 
through a self-contact mechanism [13]. For rigid scales, 
this self-contact mechanism was found to ultimately lead 
towards a kinematically dictated locked configuration. In 
the idealized system shown in Fig.j^A) where scales de¬ 
formation is assumed to be negligible with respect to the 
elastic substrate, the extent of nonlinearity in system re¬ 
sponse depends on the ratio of scale length to scale spac¬ 
ing (77 = l/d), scale embedded length (L) , scale thickness 
{D) and the elastic modulus of the substrate (F)[T3|. For 
instance, variation in 77 can be used to significantly influ¬ 
ence both bending stiffness and elastic energy absorption 
variation with curvature. Fig. Bb -C). In this context, 
the effect of friction becomes important for furthering 
the understanding of these systems since the origin of 
nonlinearity stems from scale contacts. 

In this letter we investigate the role of friction in mod- 



FIG. 1: (A) (left)A schematic of biomimetic scale sys¬ 

tem studied in the paper along with the (inset) geometri¬ 
cal dimensions (right) A 3—D printed material specimen of 
the proposed biomimetic structure (B) Non-dimensionalized 
moment-curvature response for biomimetic systems of unit 
width with various overlap ratios (77) and scales at an ini¬ 
tial angle of 10 degrees. Scale geometrical ratios assumed as 
L/D — V5,h/L — 5 and d/D > 5. The thick black curve indi¬ 
cates virgin beam with no scales with (inset) schematic of de¬ 
flection shapes plotted at nh — 0.7. (C) Non-dimensionalized 
energy density-curvature response for the same system. 


ulating the influence of biomimetic scales on a beam in 
pure bending. We assume rigid scales, linear elastic sub¬ 
strate and Coulomb type friction in order to isolate the 
role of critical geometrical and interfacial parameters in 
fostering mechanical nonlinearity as distinct from any 
traditional material sources. The nature of the nonlinear¬ 
ity becomes apparent if the system is simplified and peri¬ 
odic self-similarity of contacting scales is assumed. This 
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periodicity allows us to isolate a fundamental representa¬ 
tive volume element (RVE) for further analysing the sys- 
tern, Fig. [ 2 ]; A). A geometrical scrutiny of the RVE results 
in the following nonlinear relationship between the scale 
rotation angle and the substrate rotation, Eig. 


cos 7/^/2 — sin(6> +'0/2) = 0 (1) 

where '0 is the angular rotation of the substrate of 
RVE with respect to the instantaneous center of curva¬ 
ture of the beam, Eig. |^A) and which is directly related 
to instantaneous curvature hz of the substrate as 0 = hid. 
We also impose d <11 j2 since we are currently consider¬ 
ing only the right chirality inclination and 'ijj < a ^ the 
geometrically permissible range. Note that this relation¬ 
ship will not be satisfied if the 77 is too low, indicating 
no engagement. In general this limiting ratio is given 
by r]cr = sin Oq / Oq where Oq is the initial scale inclination 
[13]. Eor values of 77 > 77cr, the structure starts to deviate 
from the linear behavior after scale engagement according 
to Eq. The equation predicts that as the scales slide 
sufficiently from their initial inclination, they approach a 
limiting locking configuration given by Oiock and a corre¬ 
sponding limiting curvature given by i^iock where sliding 
is impossible. These terminal angles can be shown to be 
inter-related by the relationship Oiock + 0Zoc/c/2 = 7r/2 

[13]. 

In order to understand the role of friction in influenc¬ 
ing this behavior, we first look at the free body diagram 
of the RVE post-engagement, Eig. |^A). Note that the 
frictional force ffr depends on the direction of relative 
motion between scales. This can be discerned a priori 
from the kinematic relationship derived earlier, Eq. 
whose solution has two branches. Both branches meet 
at the locking configuration {Oiock^ i^iock) and depend¬ 
ing on the initial value of 0 = 6>o, one of them would 
be physically relevant. We call the branch, which corre¬ 
sponds to ^0 < Oiock ^ the low angle branch where both 
scale angle and curvature increase and the other one cor¬ 
responding to Oq > Oiock ^ the high angle branch where 
the scale angle decreases with curvature. Thus, we can 
express the friction force in terms of magnitude /o as 
ffr = fosgn{0iock — Oq) and obtain the following expres¬ 
sion from the balance of moments at the base of a single 
scale, Eig. A) 


V(/cos0-r)-/o/ sgn{0iock-0o)sm2p = Kb{0-0q) (2) 

where, r is the distance from the base to the 
point of contact between the scales, Eig. HA), Kb = 
0.66 ED'^{L/Df '^^ [T3] is the rotational stiffness mea¬ 
suring the interaction (assumed linear) of the rigid scale 
with the elastic material into which it is embedded and 
N is the inter-scale normal reaction force. The normal 
reaction force can itself be written in terms of /o and the 
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FIG. 2 : (A) Scale representative volume element (RVE) ge¬ 
ometry along with free body diagram of an individual scale 
(B) Non-dimensionalized {hie is the engagement curvature) 
frictional force for various friction coefficients. It is clear that 
friction forces approach singularity near the frictional locking 
configuration. In this simulation 77 = 6 and the initial inclina¬ 
tion Oq — 10^ (C) Phase map representation of the biomimetic 
system delineates the role of friction in mechanical behavior. 
The map shows three distinct regimes of performance linear, 
kinematically determined nonlinear and a frictionally deter¬ 
mined locking boundary (D) A schematic visual representa¬ 
tion of the locking state indexed by friction coefficient and 
77 . 

coefficient of friction fi using the Coulumb friction law 
/o ^ with the equality holding in the kinetic friction 
regime. We use the same coefficient of friction for the 
static as well as kinetic regime although in practice, the 
static friction coefficient is slightly more than its kinetic 
counterpart. Thus we arrive at the following expression 
for the magnitude of the frictional contact force (see SM 
for derivation): 

° Kb ~ ^ cos(V> + sgn{0iock - 0o)P) - r cos/3 

(3) 

where r = r/l^ ^ = tan ^ g (friction angle). Note that 
the equality is imposed as soon as sliding commences. 

The introduction of friction alters the purely kinematic 
nature of locking behavior described earlier. In this case, 
the locking mechanism is now governed by inter facial 
rather than purely kinematic variables. This can be as¬ 
certained from the frictional relationship, Eq§ which 
would yield, at a certain curvature, an infinite frictional 
force giving rise to another locking mechanism governed 
by friction. As the system approaches this locking re¬ 
gion, motion arrest begins due to escalating frictional 
forces, Eig. I^B), resulting in a limiting kinetic friction 
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FIG. 3: Non-dimensionalized (/^e is the engagement curva¬ 
ture) post engagement moment-curvature relationships for 
various values of friction coefficient clearly showing the per¬ 
ceptible difference brought about by friction in effective bend¬ 
ing stiffness of the biomimetic structure. For this simulation 
?7 = 6 and the initial inclination Oq = 10^. Scale geometrical 
ratios assumed as L/D = 15, h/L — ^ and djD > 5. 


lock which would always be lesser than the theoretical 
frictionless locking curvature. 

The overall nonlinear behavior of this frictional system 
can be unified through a phase map spanned by ^ 
which plots the dependence of scale rotation 0 with sub¬ 
strate rotation -0 for various values of r] based on Eqs|^ 
and Fig. [^C-D). In this map, the frictional locking 
configuration, which is governed by the frictional coef¬ 
ficient is represented by the kinetic frictional lock line 
in the phase space defined by + (V^/oc/c)/^ = 7 r /2 — 
sgn{0iock — where prime denotes frictional locking 
configuration. This expression quantifies the friction me¬ 
diated advance of the locking envelope alluded to earlier. 
Moreover, the nature of Coulomb friction ensures that as 
long as the initial scale angle corresponds to the low angle 
branch, any initial angle greater than the frictional lock¬ 
ing angle would cause an instantaneous post-engagement 
static frictional lock of the system. However, if the ini¬ 
tial angle ensures the higher angle branch i.e. Oq > Oiocki 
then any initial angle lower than the frictional lock an¬ 
gle will result in the static frictional lock. Thus clearly 
there exists a static friction lock band in the phase space 
of width 2 ^. Thus, a reference configuration defined by 
an initial angle and r] may cause the system to fall in 
this band and thereby rendered rigid immediately after 
engagement. However, if the configuration lies outside 


this band, the scales will move till they reach the bound¬ 
ary of this zone. Alternatively, the critical static fric¬ 
tional locking coefficient fi' for a given initial angle and 
f] can also be evaluated from the frictional locking line, 
^o + (V’Lcfe(^>M ))/2 = TT/ 2 -sgn{ 0 iock-Oo)tan~'^ iJ,. With- 
out loss of generality, taking the low angle branch, we 
can obtain the critical static frictional coefficient by solv¬ 
ing the transcendental equation /i' = cot(6>o + 6>o)). 

This shows that, for a given reference configuration, a 
frictional coefficient greater than g' would result in an 
instantaneous post-engagement locking corresponding to 
static frictional lock. Interestingly, from Fig. |^C) it be¬ 
comes clear that static frictional locking would be harder 
to achieve in the higher 77 and lower initial angles re¬ 
gion. Thus closely packed and grazing (low angle) scales, 
widely found in nature seem to be particularly resistant 
to frictional locking. 

Next, in order to understand the mechanics of this sys¬ 
tem we now derive the moment—curvature relationship of 
this structure. Assuming periodicity and Euler-Bernoulli 
beam behavior, we equate the work done due to bending 
to the total energy absorbed by the system and the fric¬ 
tional dissipation. This leads to the following variational 
energetic equation: 



Ieik^ + [l^Knie - eof+ 


1 r 

^ y /o sgn{9iock 


60) dr 


H{k- Ke) 


( 4 ) 


where M is the applied moment, k, Kq are respectively 
the current and engagement curvatures of the beam, I is 
the second moment of area of the beam cross-section and 
H{.) is the Heaviside step function for tracking engage¬ 
ment. Note that although the frictional work contains a 
singular kernel, the work integral exists. We justify lin¬ 
ear material behavior assuming relatively small strains 
in the substrate. Next, taking derivative of EqJ^ with 
respect to the current curvature yields the moment cur¬ 
vature relationship which can be evaluated numerically 
(See SM for derivation). We plot this relationship in a 
non-dimensional form for a system with scale geometry 
described in the introduction of the letter with various 
friction coefficients in the post-engagement regime to iso¬ 
late frictional contribution, Eig. This figure clearly 
shows that friction can perceptibly increase the bending 
stiffness of the structure and since it also controls the 
locking envelope, the increase can be made to occur at 
much lower curvatures. Thus friction has a dual contri¬ 
bution to the moment-curvature relationship modulating 
both the magnitude and shape of the response. 

This dual contribution of the frictional forces i.e. both 
increasing bending stiffness and lowering the locking cur¬ 
vature can have significant implications on the maximum 
frictional work during bending. To this end we define the 
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FIG. 4: Non-dimensional relative energy dissipation (RED) 
factor phase map spanned by the rj and frictional coefficient fi. 
The dissipation landscape shows maxima at the interior away 
from the boundaries indicating dissipation maximization at 
intermediate frictional coefficients. Scale geometrical ratios 
assumed as L/D = 15, h/L — ^ and djD > ^. 

relative energy dissipation (RED) factor as the ratio of 
frictional work Wfr (per unit length) to the work done on 
the system till locking Wsys (per unit length), 

RED = ^ (5) 

'^sys 

f 

where wjr = ^ /o sgn{eiock - ^o) dr, w^yg = Wfr + 

Uei with Uei = \EIk^ + 5 the elastic 
energy density of the system and the frictional locking 
curvature is = V^/oc/c/^* No'te that RED is a func¬ 
tion of 77 , friction coefficient fi and initial scale inclination 
^ 0 - However, fixing the initial angle = 10*^ for the cur¬ 
rent simulation allows us to obtain a dissipation phase 
plot which maps this non-dimensional dissipation factor 
through the system parameters 77 and /i. Fig. This 
map shows the RED factor passes through maxima with 
increasing fi and then wanes for any given 77 . Thus uni¬ 
laterally increasing fi does not necessarily improve the 
dissipation capacity of the system. This indicates that in 
the limits of low and high friction the relative contribu¬ 
tion of frictional work diminishes but through different 
mechanisms — lower overall moment contribution at low 
fi and through lowering of the locking limit for higher /i, 
corresponding to the dual nature of the frictional force. 
Finally, note that this frictional mechanism does not stem 
from or give rise to any distinction between the loading 


and unloading paths unlike classical hysteresis. 

In conclusion, we find that scale friction can advance 
the locking envelope and make high angle motion inac¬ 
cessible. Interestingly, if friction coefficient is larger than 
a critical value for a given configuration, it can instan¬ 
taneously arrest post-engagement motion, the so called 
static friction lock. Thus friction introduces a lock band 
in the system which can have useful technological im¬ 
plications. Furthermore, we find that depending on the 
sought role of the system in a particular application, such 
as stiffness gain or energy attenuation, the correspond¬ 
ing frictional property desired may be different. Thus, 
our study demonstrates that in addition to scale geom¬ 
etry and substrate properties, engineering scale surfaces 
to modulate friction can provide an important avenue for 
designing and tuning biomimetic scales systems. 
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